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Results  of  an  experimental  investigation  of  the  flow  over  a  model  roughness  are  presented.  The 
series  of  roughness  consists  of  close-packed  pyramids  in  which  both  the  height  and  the  slope  were 
systematically  varied.  The  aim  of  this  work  was  to  document  the  mean  flow  and  subsequently  gain 
insight  into  the  physical  roughness  scales  which  contribute  to  drag.  The  mean  velocity  profiles  for 
ah  nine  rough  surfaces  collapse  with  smooth-wall  results  when  presented  in  velocity-defect  form, 
supporting  the  use  of  similarity  methods.  The  results  for  the  six  steepest  surfaces  indicate  that  the 
roughness  function  A  U+  scales  almost  entirely  on  the  roughness  height  with  little  dependence  on  the 
slope  of  the  pyramids.  However,  A U+  for  the  three  surfaces  with  the  smallest  slope  does  not  scale 
satisfactorily  on  the  roughness  height,  indicating  that  these  surfaces  might  not  be  thought  of  as 
surface  “roughness”  in  a  traditional  sense  but  instead  surface  “waviness.” 

[DOI:  10.1063/1.3059630] 


I.  INTRODUCTION 

Surface  roughness  can  have  a  significant  effect  on  a  wide 
range  of  engineering  flows.  Examples  include,  but  are  not 
limited  to,  industrial  piping  systems,1  open  channel  flows,2 
turbomachines,3  marine  vehicles,4  and  aircraft.5  Because  of 
the  ubiquitous  nature  of  these  flows  in  practice,  there  has 
been  a  long  history  of  research  on  the  effect  of  roughness  on 
wall-bounded  turbulence.  Some  of  the  earliest  work  was  car¬ 
ried  out  by  Darcy6  who  investigated  the  pressure  losses  aris¬ 
ing  from  roughness  in  a  pipe.  Nikuradse7  subsequently  con¬ 
ducted  a  thorough  set  of  experiments  on  the  effects  of 
uniform  sand  roughness  on  turbulent  pipe  flow  which  signifi¬ 
cantly  furthered  the  understanding  of  roughness  effects. 
Later,  Moody,8  guided  largely  by  the  findings  of  Colebrook,9 
developed  a  diagram  to  predict  the  head  losses  in  smooth  and 
rough  pipes.  The  impact  of  the  Moody  diagram  is  difficult  to 
overstate,  as  it  has  been  a  cornerstone  in  the  held  of  hydrau¬ 
lic  engineering  for  over  60  years. 

During  the  same  period,  there  has  also  been  a  large  body 
of  research  focused  on  a  better  fundamental  understanding  of 
the  turbulence  structure  over  rough  walls.  The  classical  idea 
is  that  the  roughness  only  exerts  a  direct  influence  on  the 
turbulence  within  a  few  roughness  heights  of  the  wall,  and 
the  outer  how  is  unaffected  except  in  the  role  the  roughness 
plays  in  determining  the  outer  velocity  and  length  scales. 
This  is  often  referred  to  as  Townsend’s  hypothesis.10  The 
recent  review  of  Jimenez11  concluded  that  most  experimental 
evidence  supports  Townsend’s  hypothesis  provided  the 
roughness  is  not  too  large  compared  to  the  boundary  layer 
thickness.  Since  this  review,  the  concept  of  outer  layer  simi¬ 
larity  in  the  mean  velocity,1,12-14  Reynolds  stresses,12,15-17 
and  large  scale  turbulence  structure  ’  has  received  further 
experimental  support  for  hows  over  a  wide  range  of  three- 
dimensional  roughness. 

However,  despite  the  considerable  effort  devoted  to 
roughness  research,  many  questions  remain  unresolved.  Per¬ 


haps  the  most  important  practical  issue  is  how  to  predict  the 
frictional  drag  (or  for  internal  hows,  the  head  loss)  of  a  ge¬ 
neric  surface  based  on  measurements  of  the  roughness  topog¬ 
raphy.  At  present,  it  is  not  clear  which  roughness  length 
scales  and  parameters  best  describe  a  surface  in  a  hydraulic 
sense.  Over  the  years,  many  investigators  have  worked  on 
this  problem  (e.g.,  Refs.  19-28).  A  wide  range  of  surface 
parameters  have  been  identihed  as  correlating  with  frictional 
drag.  These  include  the  roughness  height  (e.g.,  root  mean 
square  height  ^ms,13  maximum  peak  to  trough  height  kt21 
mean  amplitude  ka ,  etc.),  slope,  density,  and  aspect 
ratio.  Higher  moments  of  the  surface  amplitude  probability 
density  function22  as  well  as  the  first  three  even  moments  of 

30 

the  wavenumber  power  spectra  of  the  surface  amplitude 
have  also  been  suggested  as  correlating  the  frictional  drag  of 
rough  surfaces.  Even  with  modest  success  of  these  correla¬ 
tions  for  a  specific  roughness  type,  it  can  be  concluded  that, 
at  present,  there  is  no  sufficiently  satisfactory  scaling  for  a 
generic,  three-dimensional  roughness.  In  fact,  the  lack  of 
meaningful  progress  toward  this  goal  led  Grigson31  to  assert 
that  the  statistics  of  the  surface  profile  alone  could  never  be 
expected  to  allow  reliable  frictional  drag  predictions,  and 
experimental  tests  of  the  surface  of  interest  would  always  be 
required.  It  is  hoped,  however,  that  investigations  which 
combine  reliable  fluid  mechanics  measurements  with  de¬ 
tailed  documentation  of  surface  topography  may  allow  fur¬ 
ther  progress  on  roughness  scaling. 

The  usual  avenue  to  studying  the  effect  of  a  given  pa¬ 
rameter  is  to  keep  ah  other  parameters  fixed  and  examine  the 
effect  of  varying  the  parameter  of  interest.  This  is  difficult 
with  surface  roughness  since  for  most  irregular  three- 
dimensional  roughness  altering  the  roughness  height,  for  ex¬ 
ample,  leads  to  concomitant  changes  in  a  host  of  other  pa¬ 
rameters.  For  this  reason,  there  are  very  few  data  available 
for  three-dimensional  roughness  in  which  the  surface  topog¬ 
raphy  is  systematically  varied.  In  the  present  experimental 
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investigation,  boundary  layer  measurements  are  made  of  the 
flow  over  a  model  roughness.  The  series  of  roughness  con¬ 
sists  of  close-packed  pyramids  in  which  both  the  height  and 
the  slope  are  varied.  The  aim  of  this  work  is  to  gain  insight 
into  how  both  roughness  height  and  slope  contribute  to  fric¬ 
tional  drag  for  three-dimensional  roughness. 


II.  BACKGROUND 


The  mean  velocity  in  the  overlap  region  of  the  inner  and 
outer  layers  for  smooth- wall-bounded  turbulent  flows  is  well 
represented  by  the  classical  log  law,32 

U+=~\n{y+)  +  B ,  (1) 

K 

where  U+  is  the  inner-normalized  mean  streamwise  velocity, 
k  is  the  von  Karman  constant  ^0.41,  y+  is  the  inner- 
normalized  distance  from  the  wall,  and  B  is  the  smooth-wall 
log-law  intercept  ^5.0.  For  flows  over  rough  walls,  there  is 
an  increased  momentum  deficit  arising  from  the  pressure 
drag  on  the  roughness  elements.  Hama  noted  that  the  pri¬ 
mary  result  of  this  is  a  downward  shift  in  the  inner- 
normalized  mean  velocity  profile  termed  the  roughness  func¬ 
tion  A  U+.  The  log  law  for  flow  over  a  rough  surface  is, 
therefore,  given  as 

U+=-ln(y+)+B-AU+.  (2) 

K 


Hama  also  observed  that  the  mean  flow  in  the  outer  layer 
was  unaffected  by  the  roughness  and  followed  the  same  pro¬ 
file  as  the  smooth  wall.  This  is  called  the  velocity-defect  law 
and  is  given  as34 


Ue-U 

UT 


(3) 


where  Ue  is  mean  streamwise  velocity  at  the  edge  of  the 
boundary  layer,  UT  is  the  friction  velocity,  and  8  is  the 
boundary  layer  thickness. 

The  roughness  function  A  U+  is  directly  related  to  the 
increase  in  skin  friction  resulting  from  the  roughness.  For 
example,  in  a  turbulent  boundary  layer,  AU+  is  a  function  of 
the  skin-friction  coefficients  Cj  for  the  smooth  and  rough 
walls  at  the  same  displacement  thickness  Reynolds  number 
Re^  given  as33 


The  roughness  function  depends  on  both  the  nature  of 
the  roughness  and  the  Reynolds  number  of  the  flow  such  that 

A  U+=f(k+).  (5) 


Once  the  functional  dependence  expressed  in  Eq.  (5)  is 
known  for  a  given  rough  surface,  the  frictional  drag  for  an 
arbitrary  body  covered  with  that  roughness  can  be  predicted 
at  any  Reynolds  number  using  a  computational  boundary 
layer  code  or  a  similarity  law  analysis. 

One  difficulty,  however,  is  that  the  roughness  function  is 
not  universal  among  roughness  types.  This  is  clearly  illus- 
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•  Floned  Surface  -  BL  [12] 

o  Floned  Surface  -  Pipe  [1 3] 

♦  Commerical  Steel  Pipe  [1] 

.  Uniform  Sand  [7,  36] 

-  Colebrook-type  RF  [9] 

- Fully-Rough  Asymptote 
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FIG.  1 .  Comparison  of  the  roughness  functions  measured  in  some  previous 
studies. 


trated  in  Fig.  1.  The  results  presented  are  for  honed 
surfaces,  ’  commercial  steel  pipe,  and  uniform  sand.  ’ 
Also  shown  is  the  Colebrook-type9  roughness  function  for 
“naturally  occurring”  surfaces,  upon  which  the  Moody  dia¬ 
gram  is  based.  The  most  obvious  disparity  in  these  roughness 
functions  is  their  behavior  in  the  transitionally  rough  regime. 
The  critical  roughness  Reynolds  number  k+  for  both  the  on¬ 
set  of  roughness  effects  and  the  approach  to  the  fully  rough 
behavior  depends  on  the  roughness  type.  A  more  fundamen¬ 
tal  issue  that  has  plagued  researchers  is  how  to  relate  the 
roughness  topography  to  k ,  the  hydraulic  roughness  length 
scale.  That  is,  what  roughness  length  scale  will  cause  col¬ 
lapse  of  the  roughness  functions  for  a  range  of  surfaces  in  the 
fully  rough  flow  regime?  Typically  the  length  scale  which  is 
adopted  is  the  equivalent  sand  roughness  height  ks,  but  its 
relation  to  the  topography  of  a  generic  surface  is  unclear.  ks 
(Ref.  36)  is  defined  as  the  size  of  uniform  sand  in  Ni- 
kuradse’s  experiments  that  gives  the  same  A  U+  as  the  surface 
of  interest  in  the  fully  rough  flow  regime.  This  relationship 
can  be  expressed  as36 

AU+  =  — ln(&*)  +  7?  -  8.5.  (6) 

K 

If,  for  example,  one  examines  the  surfaces  in  Fig.  1,  ks  is 
simply  the  size  of  the  sand  grains  for  the  uniform  sand  sur¬ 
face,  while  ks~3krms  for  the  honed  surfaces12,13  and  ks 
~  1.6/:^  for  the  commercial  steel  pipe.1  Clearly  the  rough¬ 
ness  height  alone  (k^  in  this  example)  is  not  sufficient  to 
specify  ks  for  a  given  surface.  In  the  present  investigation,  it 
was  hypothesized  that  the  slope  of  the  roughness  elements 
should  also  be  important  in  determining  ks.  This  is  plausible 
since  separation  of  the  flow  off  of  individual  roughness  ele¬ 
ments  should  depend  on  the  slope.  Previous  work  by 
Musker  supports  this  view.  The  use  of  a  model,  close- 
packed  pyramid  surface  allows  the  roughness  height  and 
slope  to  both  be  varied  while  other  statistics  such  as  the 
skewness  and  the  kurtosis  of  the  roughness  amplitude  remain 
fixed.  In  this  way,  it  is  hoped  that  the  role  of  the  roughness 
height  and  slope  in  determining  ks  can  be  more  clearly  elu- 
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FIG.  2.  Schematic  of  the  flat  plate  test  fixture. 

FIG.  3.  Schematic  of  the  roughness  geometry  and  arrangement  illustrating 
the  roughness  height  kt  and  slope  angle  a. 


cidated.  This  work  is  part  of  a  larger  effort  to  collect  results 
for  three-dimensional  rough  surfaces  in  an  effort  to  deter¬ 
mine  the  roughness  scales  which  contribute  to  frictional 
drag. 

III.  EXPERIMENTAL  FACILITIES  AND  METHOD 

The  experiments  were  conducted  in  the  high-speed  water 
tunnel  facility  at  the  United  States  Naval  Academy  Hydro¬ 
mechanics  Laboratory.  The  test  section  has  a  cross  section  of 
40  X  40  cm2  and  is  1.8  m  in  length.  A  velocity  range  of  0—9 
m/s  can  be  produced  in  the  test  section.  In  the  present  ex¬ 
periments,  four  freestream  velocities  were  tested  ranging 
from  1  to  7  m/s.  Flow  management  devices  in  the  facility 
include  turning  vanes  in  the  tunnel  corners  and  a  honeycomb 
flow  straightener  in  the  settling  chamber.  The  honeycomb 
has  19  mm  cells  that  are  150  mm  in  length.  The  area  ratio 
between  the  settling  chamber  and  the  test  section  is  20: 1 .  The 
resulting  freestream  turbulence  intensity  in  the  test  section  is 
—0.5%. 

The  test  surfaces  were  mounted  into  a  splitter-plate  test 
fixture.  The  test  fixture  was  installed  horizontally,  at  mid¬ 
depth  in  the  tunnel.  The  first  0.20  m  of  the  test  fixture  was 
covered  with  No.  36-grit  sandpaper  to  fix  the  location  of 
transition  and  thicken  the  resulting  turbulent  boundary  layer. 
A  schematic  of  the  test  plate  fixture  is  shown  in  Fig.  2.  Fur¬ 
ther  details  of  the  experimental  facility  can  be  found  in  Ref. 
12.  All  the  measurements  reported  here  were  obtained  at  v 
=  1.35  m  downstream  of  the  leading  edge.  Velocity  profiles 
taken  upstream  of  the  measurement  location  indicated  that 
self- similarity  in  velocity-defect  profiles  was  achieved  for  all 
the  test  surfaces  for  x>0.80  m.  The  upper  removable  wall 
of  the  tunnel  is  adjustable  to  account  for  boundary  layer 
growth.  In  the  present  work,  this  wall  was  set  to  produce  a 
nearly  zero  pressure  gradient  boundary  layer.  The  strength  of 
the  pressure  gradient  can  be  quantified  by  the  acceleration 
parameter  K  given  as 


square,  right  pyramids.  The  four  lateral  edges  of  the  pyra¬ 
mids  were  oriented  at  0°,  180°,  and  ±90°  to  the  freestream 
flow  (Fig.  3).  Three  pyramid  heights  kt  were  tested  with  kt 
~0.30,  0.45,  and  0.60  mm.  The  slope  angle  of  the  lateral 
edges  with  the  horizontal,  a,  was  also  varied.  Slope  angles  of 
a~ll°,  22°,  and  45°  were  tested.  The  combination  of  the 
three  pyramid  heights  and  three  slope  angles  produced  the 
nine  rough  surfaces  that  were  tested  in  the  present  study  (Fig. 
4).  The  series  of  test  surfaces  were  produced  using  three 
different  custom  carbide  engraving  cutters,  one  for  each  de¬ 
sired  slope  angle.  The  surfaces  were  machined  using  a  CNC 
Haas  VF-11,  three-axis  vertical  milling  machine,  and  ESPRIT 
CAD/CAM  software.  The  machining  was  carried  out  at  ±45° 
to  the  direction  of  the  flow.  The  tool  path  pitch  and  depth  of 
cut  were  controlled  to  produce  the  desired  surface  texture. 

In  order  to  ensure  the  fidelity  of  the  fabrication,  three- 
dimensional  topographical  profiles  were  made  of  the  rough 
surfaces  using  a  Micro  Photonics  Nanavea  ST300  white  light 
chromatic  aberration  surface  profilometer.  The  vertical  accu¬ 
racy  of  the  system  is  0.3  jum  with  lateral  resolution  of 
6  fim.  The  data  were  digitized  at  increments  of  25  jmm  in 
the  lateral  directions,  and  the  sampling  area  was  sufficient  to 
capture  a  few  roughness  elements.  An  example  of  a  surface 
map  for  the  /q~0.45  mm,  11°  roughness  is  presented  in 
Fig.  5.  The  topographical  profiles  indicate  that  the  roughness 
height  was  within  ±5%  of  the  desired  height  for  all  the  sur¬ 
faces,  while  the  slope  angle  was  within  ±1%.  The  experi¬ 
mental  test  conditions  along  with  the  roughness  designation 
used  are  given  in  Table  I. 

Velocity  measurements  were  obtained  using  a  TSI 
FSA3500  two-component  laser  Doppler  velocimeter  (LDV). 
The  LDV  system  utilized  a  four  beam  fiber  optic  probe  and 
was  operated  in  backscatter  mode.  The  probe  included  a  cus- 


v  dUe 

K=~; - £ 

u  :  dx 


(7) 


In  the  present  experiments,  K<  1  X  10-8  in  all  cases,  which 
is  more  than  an  order  of  magnitude  below  the  value  where 
pressure  gradient  effects  are  likely  to  influence  the  mean 
flow.37 

The  removable  test  specimens  were  fabricated  from  12 
mm  thick  cast  acrylic  sheet  of  350  mm  in  width  and  1.32  m 
in  length.  The  rough  surfaces  consisted  of  close-packed, 


a  =  45°  a  =  22° 


a-  11° 


k,  =  0.30  mm  A 


kt  =  0.45  mm 


kt  =  0.60  mm  A 


FIG.  4.  Matrix  of  roughness  test  surfaces.  The  roughness  is  not  drawn  actual 
size  but  is  to  scale. 
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0  2  4  6  8  10  mm  mm 


mm 


FIG.  5.  (Color  online)  Plan  view  of  the  roughness  topography  for  the  a 
=  22°,  kt= 450  /mm  specimen  measured  with  a  profilometer  (overall  uncer¬ 
tainty  in  roughness  elevation:  ±3  X  10-4  mm). 


tom  beam  displacer  to  shift  one  of  the  four  beams.  The  result 
was  three  coplanar  beams  aligned  parallel  to  the  wall  along 
with  a  fourth  beam  in  the  vertical  plane  which  allowed  mea¬ 
surements  near  the  wall  without  having  to  tilt  or  rotate  the 
probe.  The  system  also  employed  2.6:1  beam  expansion  op¬ 
tics  at  the  exit  of  the  probe  to  reduce  measurement  volume 
size.  The  resulting  probe  volume  diameter  d  was  45  jum, 
with  a  length  /  of  340  jum.  The  measurement  volume  diam¬ 
eter,  therefore,  ranged  from  1.9<d+<16,  while  its  length 
was  14</+<120  when  expressed  in  viscous  length  scales. 
The  range  of  d+  values  in  the  present  experiment  is  similar  to 
the  smooth- wall  turbulent  boundary  layer  study  of  DeGraaff 
and  Eaton.38 

In  the  present  investigation,  20  000  random  velocity  re¬ 
alizations  were  obtained  at  each  location  in  the  boundary 
layer,  and  the  data  were  collected  in  coincidence  mode.  A 
Velmex  three-axis  traverse  unit  allowed  the  position  of  the 
LDV  probe  to  be  maintained  to  ±5  jmm  in  all  directions. 
Velocity  profiles  consisted  of  approximately  42  sampling  lo¬ 
cations  within  the  boundary  layer.  The  flow  was  seeded  with 
3  jum  diameter  alumina  particles.  The  seed  volume  was  con¬ 
trolled  to  achieve  acceptable  data  rates  while  maintaining  a 
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low  burst  density  signal. 

The  friction  velocity  was  determined  using  two  methods. 
The  first  was  the  total  stress  method.40  This  method  assumes 
that  the  plateau  in  the  total  shear  stress  in  the  inner  layer  is 
equivalent  to  the  shear  stress  at  the  wall.  The  total  shear 
stress  was  calculated  by  adding  the  viscous  and  turbulent 
stress  components.  The  second  method  employed  for  finding 
UT  was  the  modified  Clauser  chart  method.41  Details  of  both 
methods  used  to  find  the  wall  shear  stress  can  be  found  in 
Ref.  17.  The  normalization  presented  throughout  the  remain¬ 
der  of  this  paper  uses  the  friction  velocity  obtained  using  the 
modified  Clauser  chart  method.  It  should  be  noted  that  the 
friction  velocity  values  found  using  modified  Clauser  chart 
and  total  stress  methods  agreed  within  <3.0%  for  all  cases 
tested  with  a  mean  absolute  difference  of  1.3%. 


IV.  UNCERTAINTY  ESTIMATES 

Estimates  of  the  overall  uncertainty  in  the  measured 
quantities  were  made  by  combining  their  precision  and  bias 
uncertainties  using  the  methodology  specified  by  Moffat.42 
The  95%  precision  confidence  limits  for  a  given  quantity 
were  obtained  by  multiplying  its  standard  error  by  the  two- 
tailed  t  value  given  by  Coleman  and  Steele.43  Standard  error 
estimates  were  determined  from  the  variability  observed  in 
repeated  velocity  profiles  taken  on  a  given  test  surface.  Bias 
estimates  were  combined  with  the  precision  uncertainties  to 
calculate  the  overall  uncertainties  for  the  measured  quanti¬ 
ties.  It  should  be  noted  that  the  LDV  data  were  corrected  for 
velocity  bias  by  employing  burst  transit  time  weighting.44 
Fringe  bias  was  considered  insignificant,  as  the  beams  were 
shifted  well  above  a  burst  frequency  representative  of  twice 
the  freestream  velocity.45  The  resulting  overall  uncertainty  in 
the  mean  velocity  was  ±1.5%.  The  overall  uncertainty  in  UT 
obtained  from  the  modified  Clauser  chart  method  was  ±4%, 
while  the  total  stress  method  yielded  an  uncertainty  of  ±6%. 

V.  RESULTS  AND  DISCUSSION 

The  presentation  of  the  results  and  discussion  will  be 
organized  as  follows.  First,  the  mean  flow  profiles  over  the 
systematically  varied  roughness  will  be  discussed  in  terms  of 
both  inner  and  outer  scalings.  Next,  the  roughness  function 
scaling  with  regards  to  roughness  height  and  slope  will  be 
considered. 

A.  Mean  velocity  profiles 

The  mean  velocity  profiles  of  the  rough  walls  are  pre¬ 
sented  using  inner  scaling  in  Fig.  6.  For  the  sake  of  clarity 
and  brevity,  most  of  the  results  presented  in  this  section  are 
for  the  ^=450  mm  cases  only,  although  the  trends  for  the 
other  roughness  heights  are  similar.  Also  shown  for  compari- 
son  are  the  smooth-wall  results  of  DeGraaff  and  Eaton  at 
Re6>=5100.  It  can  be  seen  that  all  of  the  roughness  slope 
angles  exhibit  an  increasing  roughness  function  (Af/+)  with 
increasing  Reynolds  number.  One  of  the  most  striking  obser¬ 
vations  is  that  the  mean  velocity  profiles  for  the  two  steepest 
slope  angles  (a=22°  and  a- 45°)  are  nearly  identical  [Figs. 
6(a)  and  6(b)],  suggesting  that  the  slope  of  the  present  rough¬ 
ness  elements  has  little  influence  on  the  resulting  drag.  How¬ 
ever,  the  surface  with  the  shallowest  slope  (a^ll0)  clearly 
shows  a  smaller  A  U+  and  hence  has  less  drag  than  the  other 
surfaces.  The  roughness  slope,  therefore,  can  affect  the  drag, 
although  slope  dependence  is  only  apparent  for  shallower 
angles.  This  will  be  discussed  in  greater  detail  in  the  rough¬ 
ness  function  results. 

The  outer- scaled  mean  velocity  profiles  are  presented  in 
velocity-defect  form  in  Fig.  7.  Again,  the  smooth- wall  data 
from  DeGraaff  and  Eaton  at  Re^SlOO  are  also  shown.  All 
the  mean  velocity  profiles  exhibit  agreement  well  within  ex¬ 
perimental  uncertainty  in  the  log  law  and  outer  regions  of  the 
boundary  layer.  The  Clauser-Rotta46  length  scale  (A 
=  U+e8)  has  been  offered  as  an  alternative  outer  length  scale 
for  the  mean  velocity  profile.  It  has  the  advantage  of  being 
based  on  an  integral  length  parameter  (8)  rather  than  S, 
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TABLE  I.  Experimental  test  conditions. 


Designation 

K 

(gm) 

a 

(°) 

V. 

(m  s'1) 

Re„ 

e 

(gm) 

UT  Clauser  chart 
(m  s-1) 

UT  Total  stress 
(m  s'1) 

S 

(mm) 

K 

A  U+ 

45_300_1 

300 

45 

1.01 

3  430 

200 

0.0432 

0.0432 

28.1 

1  180 

12.6 

1.7 

45_300_3 

300 

45 

3.01 

11450 

80 

0.137 

0.138 

30.3 

3  960 

39.2 

6.2 

45_300_5 

300 

45 

4.99 

20  850 

10 

0.231 

0.233 

33.3 

7  340 

66.1 

7.8 

45_300_7 

300 

45 

7.01 

30  610 

70 

0.328 

0.328 

34.1 

10  850 

95.5 

9.1 

45_450_1 

450 

45 

1.01 

3  510 

280 

0.0458 

0.0455 

28.5 

1  250 

19.8 

3.3 

45_450_3 

450 

45 

3.01 

11  570 

60 

0.142 

0.145 

29.9 

4  130 

62.2 

7.0 

45_450_5 

450 

45 

5.00 

20  060 

150 

0.240 

0.242 

31.2 

7  140 

103.0 

00 

bo 

45_450_7 

450 

45 

7.01 

28  580 

200 

0.339 

0.334 

31.7 

10  220 

145.2 

9.8 

45_600_1 

600 

45 

1.00 

3  960 

220 

0.0498 

0.0496 

30.1 

1  470 

29.3 

5.4 

45_600_3 

600 

45 

3.01 

12  320 

50 

0.159 

0.159 

31.0 

4  740 

91.7 

9.3 

45_600_5 

600 

45 

5.00 

21  760 

150 

0.254 

0.253 

32.2 

7  930 

147.7 

10.3 

45_600_7 

600 

45 

7.01 

30  100 

200 

0.354 

0.355 

31.8 

10  920 

206.0 

10.9 

22_300_1 

300 

22 

1.01 

3  430 

200 

0.0448 

0.0447 

28.6 

1  240 

13.0 

2.6 

22_300_3 

300 

22 

3.01 

11  280 

100 

0.138 

0.135 

30.2 

4  020 

39.9 

6.4 

22_300_5 

300 

22 

5.01 

21  020 

170 

0.234 

0.227 

32.2 

7  300 

68.0 

8.3 

22_300_7 

300 

22 

7.02 

28  610 

180 

0.326 

0.318 

31.6 

9  990 

94.9 

9.0 

22_450_1 

450 

22 

1.01 

3  460 

200 

0.0467 

0.0461 

27.3 

1  240 

20.4 

3.5 

22_450_3 

450 

22 

3.01 

10  400 

40 

0.146 

0.149 

28.0 

3  900 

62.7 

7.0 

22_450_5 

450 

22 

5.00 

17  730 

10 

0.246 

0.248 

29.0 

6  850 

106.4 

8.3 

22_450_7 

450 

22 

7.01 

25  510 

130 

0.346 

0.342 

29.2 

9  800 

151.5 

9.5 

22_600_1 

600 

22 

1.01 

3  770 

240 

0.0493 

0.0493 

28.7 

1  370 

28.7 

5.1 

22_600_3 

600 

22 

3.00 

11  450 

250 

0.155 

0.151 

29.4 

4  290 

87.5 

8.7 

22_600_5 

600 

22 

5.00 

19  780 

260 

0.257 

0.250 

30.6 

7  500 

147.1 

9.9 

22_600_7 

600 

22 

7.00 

28  030 

300 

0.363 

0.357 

30.5 

10  740 

211.3 

10.8 

11_300_1 

300 

11 

1.01 

3  070 

100 

0.0429 

0.0423 

26.3 

1  090 

12.5 

1.1 

11_300_3 

300 

11 

3.00 

9  180 

70 

0.127 

0.127 

26.2 

3  140 

35.9 

3.8 

11_300_5 

300 

11 

5.00 

15  820 

60 

0.214 

0.215 

26.9 

5  490 

61.2 

5.3 

11_300_7 

300 

11 

7.02 

21  590 

100 

0.303 

0.298 

26.5 

7  800 

88.2 

6.1 

11_450_1 

450 

11 

1.00 

2  910 

130 

0.0435 

0.0445 

25.8 

1  080 

18.8 

1.5 

11_450_3 

450 

11 

3.01 

8  640 

120 

0.127 

0.129 

25.2 

3  030 

54.2 

3.9 

11_450_5 

450 

11 

5.01 

14  870 

270 

0.217 

0.221 

26.5 

5  450 

92.8 

5.3 

11_450_7 

450 

11 

7.02 

19  980 

320 

0.307 

0.305 

25.2 

7  440 

132.7 

6.2 

11_600_1 

600 

11 

1.00 

3  480 

340 

0.0429 

0.0439 

28.3 

1  180 

25.0 

2.2 

11_600_3 

600 

11 

3.00 

10  710 

400 

0.125 

0.126 

29.5 

3  550 

72.1 

4.2 

11_600_5 

600 

11 

5.01 

18  620 

400 

0.208 

0.207 

31.2 

6  290 

120.9 

5.2 

11_600_7 

600 

11 

7.01 

25  670 

480 

0.296 

0.288 

31.0 

8  940 

173.0 

6.2 

whose  definition  is  somewhat  arbitrary.  It  also  can  account 
for  Reynolds  number  dependence  through  U\.  The  outer- 
scaled  mean  velocity  profiles  are  presented  in  velocity-defect 
form  using  the  Clauser-Rotta  scaling  in  Fig.  8.  The  results 
show  excellent  collapse  throughout  the  log  law  and  outer 
regions  of  the  boundary  layer  for  both  the  smooth  and  rough 
walls.  These  mean  flow  results  support  Townsend’s 
hypothesis10  that  the  surface  condition  has  no  direct  effect  on 
the  outer  flow  but  simply  plays  a  role  in  setting  the  outer 
length  and  velocity  scales.  This  is  in  agreement  with  earlier 
work  by  the  authors12,14  for  range  of  three-dimensional 
roughness  of  various  heights  and  types.  Recent  results47 
show  that  this  mean  flow  similarity  is  very  robust  and  is 
maintained  for  flows  over  three-dimensional  roughness  in 
which  the  roughness  height  is  a  significant  fraction  of  the 
boundary  layer  thickness  (relative  roughness  as  large  as 


kt/  S~  0.2).  Collapse  of  mean  flow  profiles  for  rough  and 
smooth  walls  in  velocity-defect  form  is  important  in  accurate 
prediction  of  drag,  since  mean  flow  similarity  in  outer  vari¬ 
ables  is  assumed  both  when  relating  laboratory  results  to  full 

scale  using  boundary  layer  similarity  laws  as  well  as  when 

35 

using  wall  function  models  for  numerical  computations. 

B.  Roughness  function  AU+ 

The  roughness  function  results  for  the  two  steepest 
roughness  slope  angles  (a=22°  and  a- 45°)  are  shown  in 
Fig.  9  along  with  the  results  for  Nikuradse’s7,36  uniform 
sand.  For  the  present  results,  the  roughness  height  k  was 
initially  taken  to  be  the  peak-to-trough  height  kt.  Figure  9 
shows  reasonable  agreement  between  the  roughness  func¬ 
tions  observed  on  the  a -22°  and  a= 45°  surfaces.  Any  dif- 
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FIG.  6.  Mean  velocity  profiles  in  inner  variables  (overall  uncertainty  in 
U+:  ±4%). 


FIG.  7.  Mean  velocity  profiles  in  outer  variables  using  classical  scaling 
(overall  uncertainty  in  U+e-U+\  ±5%) 
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FIG.  8.  Mean  velocity  profiles  in  outer  variables  using  the  Clauser-Rotta 
length  scale  (overall  uncertainty  in  U+e-U+\  ±5%). 


ferences  that  are  observed  do  not  indicate  any  clear  influence 
of  the  slope  angle.  It  appears,  therefore,  if  the  roughness 
slope  is  steep  enough,  the  roughness  function  scales  entirely 
on  k ,  and  the  roughness  slope  is  irrelevant.  These  roughness 
slope  angles  combined  with  the  sharp  edges  of  the  present 
roughness  most  likely  cause  total  flow  separation  in  the  lee 
of  the  individual  roughness  elements  regardless  if  the  slope 
angle  is  22°  or  45°,  leading  to  large  form  drag. 

The  results  at  the  higher  Reynolds  numbers  were  used  to 
evaluate  the  equivalent  sand  roughness  height  ks  using  Eq. 
(6).  For  the  a =22°  and  a- 45°  surfaces,  this  gives  ks 
~  1.5 kt.  The  results  for  the  steeper  roughness  slopes  are  plot¬ 
ted  in  Fig.  10  using  this  scaling.  It  can  be  seen  that  the 
roughness  function  for  these  surfaces  displays  a  similar  in¬ 
flectional  behavior  in  the  transitionally  rough  regime  as  the 
uniform  sand.  However,  it  appears  that  for  the  present  rough¬ 
ness,  the  inflectional  behavior  is  even  more  pronounced  than 


FIG.  9.  Roughness  function  results  for  the  a=22°  and  cr=45°  specimens 
using  k=kt  (overall  uncertainty  in  A U+:  ±  10%). 
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k+ 


FIG.  10.  Roughness  function  results  for  the  a =22°  and  a =45°  specimens 
using  k=1.5kt  (overall  uncertainty  in  A U+:  ±  10%). 

is  observed  for  the  uniform  sand.  More  low  roughness  Rey¬ 
nolds  number  data  would  be  useful  to  confirm  this  observa¬ 
tion. 

The  roughness  function  results  for  the  shallowest  rough¬ 
ness  slope  angle  (a,=  ll°)  are  shown  in  Fig.  11  along  with 
the  data  for  uniform  sand  of  Nikuradse.7,36  Note  that  for  the 
present  results,  the  roughness  height  k  was  taken  to  be  the 
peak-to-trough  height  kt.  The  results  clearly  indicate  that  the 
roughness  function  for  the  a- 11°  surfaces  does  not  scale  on 
k.  For  example,  examination  of  the  higher  Reynolds  number 
results  reveals  that  doubling  the  roughness  height  has  little,  if 
any,  effect  on  A  U+.  This  seems  to  indicate  that  the  separation 
in  the  lee  of  the  shallow  slope  roughness  is  not  complete  as 
is  the  case  of  the  steeper  slopes.  That  is  not  to  say  that  there 
is  no  separation  behind  the  elements.  Both  the  behavior  of 
the  roughness  function  and  the  sharp-edged  nature  of  the 
roughness  elements  would  point  to  at  least  mild  separation. 
The  fact  that  there  is  a  positive  A  U+  for  these  surfaces  indi¬ 
cates  that  their  skin  friction  is  above  smooth- wall  values  and 
there  is  some  form  drag  on  the  roughness  elements.  Attempts 
were  made  to  use  alternative  length  scales  to  collapse  the 


FIG.  11.  Roughness  function  results  for  the  a=  11°  specimens  using  k=kt 
(overall  uncertainty  in  A  U+:  ±  10%). 


roughness  function  for  these  surfaces.  This  met  with  little 
success  because  of  the  geometric  similarity  of  the  surfaces. 
If,  for  example,  the  roughness  wavelength  is  used,  the  result 
is  effectively  the  same  as  using  k,  since  the  wavelength  of 
the  roughness  varies  in  the  same  fashion  as  the  roughness 
height  since  the  slope  of  the  surfaces  is  the  same. 

Rough  surfaces  in  which  the  roughness  function  does  not 
scale  on  k  have  been  reported  on  extensively  in  the  literature 
(e.g.,  Refs.  41  and  49-51).  However,  most,  if  not  all  of  this 
work,  has  focused  on  the  flow  over  two-dimensional  rough¬ 
ness,  specifically,  transverse  bars  in  which  the  spacing  be¬ 
tween  bars  is  nearly  the  same  as  the  bar  height.  Perry  et  al.52 
first  showed  that  this  type  of  roughness  displayed  an  anoma¬ 
lous  roughness  function  scaling  and  termed  it  “d-type” 
roughness.  This  is  in  contrast  to  most  roughness  types  in 
which  A  U+  scales  in  some  fashion  on  k ,  commonly  called 
&-type  roughness.  In  the  case  of  two-dimensional  bars,  the 
reason  for  the  anomalous  roughness  function  scaling  stems 
from  the  fact  that  the  flow  skims  over  the  cavity  between  the 
bars  leading  to  the  formation  of  stable  vortical  flow  cells 
between  elements.  To  the  authors’  knowledge  the  inability  of 
the  roughness  function  to  scale  on  the  roughness  height  has 
not  been  previously  documented  in  the  literature  for  any  type 
of  three-dimensional  roughness.  The  mechanism  for  the  lack 
of  k  scaling  in  the  present  three-dimensional  roughness  is 
obviously  much  different  than  for  the  case  of  two- 
dimensional  bars.  For  the  present  roughness,  there  appears  to 
be  a  critical  slope  below  which  the  roughness  no  longer  be¬ 
haves  like  “roughness”  in  a  traditional  sense  but  instead  acts 
as  surface  “waviness.” 

Nakato  et  al.  noted  analogous  behavior  for  two- 
dimensional  sinusoidal  wave  roughness.  They  observed  that 
the  roughness  function  behaved  like  that  for  uniform  sand 
provided  the  slope  was  greater  than  —6°.  In  cases  where  the 
slope  was  <6°,  the  roughness  function  did  not  follow  the 
uniform  sand  roughness  function,  and  they  referred  to  these 
surfaces  as  “wavy”  walls.  They  did  not  note  any  lack  of 
scaling  on  k  for  these  surfaces.  However,  this  may  have  been 
simply  because  they  did  not  systematically  vary  k  in  the 
experiments  while  holding  the  slope  fixed.  This  is  an  impor¬ 
tant  point  because,  unlike  Nikuradse’ s  exhaustive  study  in 
which  both  the  Reynolds  number  and  roughness  height  were 
varied,  many  studies  only  vary  the  Reynolds  number  while 
holding  the  roughness  height  fixed.  If  one  examines  the  re¬ 
sults  presented  in  Fig.  11  for  any  single  roughness  height,  it 
would  likely  be  concluded  that  the  roughness  function  does 
scale  on  k  since  an  increase  in  the  roughness  Reynolds  num¬ 
ber  leads  to  an  increasing  roughness  function.  However, 
when  the  results  are  taken  in  their  entirety,  it  is  obvious  this 
is  not  the  case.  There  may,  therefore,  be  more  surfaces  that 
have  been  identified  as  k-type  based  on  results  from  a  single 
roughness  height  that  are  not. 

The  present  results,  taken  in  light  of  the  recent  study  of 
Napoli  et  al.,54  do  provide  some  physical  insight  into  this 
behavior.  Napoli  et  al.  carried  out  a  numerical  investigation 
of  the  flow  over  irregularly  distributed  two-dimensional 
roughness.  The  wall  consisted  of  corrugated  roughness  con¬ 
structed  by  the  superposition  of  sinusoids  of  random  ampli- 
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FIG.  12.  Roughness  function  vs  ES  (overall  uncertainty  in  A U+:  ±  10%). 


tude.  In  this  work,  they  developed  a  roughness  parameter 
termed  the  effective  slope  (ES),  defined  as  follows: 


(8) 


where  L  is  the  sampling  length,  r  the  roughness  amplitude, 
and  v  the  streamwise  direction.  In  their  study,  roughness  with 
effective  slopes  in  the  range  of  0.042  <ES<  0.760  were 
tested,  and  it  was  concluded  that  the  roughness  function  is 
strongly  dependent  on  ES.  Specially,  they  concluded  that 
A U+  scales  linearly  on  ES  for  ES  <  0.15.  For  larger  values  of 
ES,  A U+  scales  on  ES  in  a  nonlinear  fashion,  reaching  a 
maximum  at  ES  — 0.55  and  declining  slightly  for  ES>0.55. 
It  should  be  noted  that  all  of  their  test  cases  were  carried  out 
at  the  same  friction  Reynolds  number.  The  change  in  the 
scaling  which  they  contend  occurs  at  ES  — 0.15  corresponds 
roughly  to  the  point  where  frictional  and  form  drag  on  the 
roughness  elements  are  equal. 

In  the  present  work,  the  range  of  effective  slopes  tested 
was  0.19<ES<  1.0.  The  results  of  Napoli  et  al. 54  are  pre¬ 
sented  in  Fig.  12  along  with  the  present  results  from  the 
highest  Reynolds  number  cases.  The  trends  observed  in  the 
results  are  similar.  The  present  results  show  that  for  the 
steepest  roughness  cases  (ES  =  0.40  or  1.0),  the  roughness 
function  does  not  scale  on  ES.  The  variation  in  A U+  ob¬ 
served  is  due  to  differences  in  k,  and  the  roughness  function 
scales  entirely  on  the  roughness  Reynolds  number  [Eq.  (5)]. 
This  corresponds  to  the  roughness  flow  regime  in  which 
form  drag  on  the  roughness  elements  is  much  larger  than  the 
frictional  drag.  For  the  smallest  slope  tested  (ES  =  0.19),  A U+ 
is  reduced  compared  to  higher  ES  cases  and  A  U+  does  not 
scale  on  k.  This  corresponds  to  the  waviness  flow  regime  in 
which  ES  is  an  important  parameter  in  scaling  the  roughness 
function.  In  this  regime,  the  frictional  drag  is  non-negligible. 
Based  on  the  present  results  and  those  of  Napoli  et  al.,54  it  is 
proposed  that  a  regime  transition  occurs  at  ES  — 0.35.  For 
ES>— 0.35  (the  roughness  flow  regime),  the  roughness 
function  is  independent  of  the  effective  slope  and  scales  en¬ 
tirely  on  the  roughness  Reynolds  number.  For  ES<  —0.35 
(the  roughness  flow  regime),  the  roughness  function  is 


strongly  dependent  on  the  effective  slope,  while  the  rough¬ 
ness  height  has  little  influence.  In  light  of  the  results  of 
Napoli  et  al.54  this  transition  appears  to  correspond  physi¬ 
cally  to  the  slope  at  which  pressure  drag  on  the  roughness 
elements  overwhelms  the  frictional  drag.  However,  a  detailed 
survey  of  the  flow  structure  very  near  the  roughness  elements 
is  needed  to  better  understand  the  flow  physics  which  give 
rise  to  the  transition  from  the  roughness  to  the  waviness  re¬ 
gime. 


VI.  CONCLUSION 

Results  of  an  experimental  investigation  of  the  flow  over 
a  model  roughness  in  which  both  the  height  and  the  slope 
were  systematically  varied  have  been  presented.  The  mean 
velocity  profiles  for  all  the  rough  surfaces  collapse  with 
smooth-wall  results  when  presented  in  velocity-defect  form. 
The  slope  of  the  roughness  appears  to  be  an  important  pa¬ 
rameter  in  the  prediction  of  drag  for  roughness  with  shallow 
angles;  however,  its  importance  diminishes  with  increasing 
slope  until  the  point  where  the  roughness  height  is  the  domi¬ 
nant  scaling  parameter.  The  results  for  the  steepest  surfaces 
indicate  that  the  roughness  function  A  U+  scales  almost  en¬ 
tirely  on  the  roughness  height  with  little  dependence  on  the 
slope  of  the  pyramids.  However,  A U+  for  the  three  surfaces 
with  the  smallest  slope  does  not  scale  satisfactorily  on  the 
roughness  height,  indicating  that  these  surfaces  might  not  be 
thought  of  as  surface  roughness  in  a  traditional  sense  but 
instead  surface  waviness.  There  appears  to  be  a  critical  ef¬ 
fective  slope,  (ES  — 0.35)  that  delineates  the  flow  regimes. 
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